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Abstract 

Ultra -compact  mode-locked  CO2  lasers  and  mode- 
locked  Ndrfiber  lasers  are  being  developed  for 
Theater  Missile  Defense  (TMD)  and  National  Missile 
Defense  (NMD)  applications.  These  lasers  produce 
a  coherent  train  of  narrow  pulses  which  may  be  used 
in  a  coherent  laser  radar  (ladar)  to  form  an  image  of 
a  possible  target.  The  distance  to  the  target  is 
irrelevant  as  long  as  there  is  enough  received  energy. 
These  images  are  known  as  "range-Doppler,"  ’Tange- 
resolved  Doppler,”  or  "inverse  synthetic  aperture 
radar  (ISAR)"  images.  The  down-range  image 
resolution  is  given  by  the  micro-pulse  width 
multiplied  by  c/2,  and  the  image  Doppler  resolution 
is  given  by  the  reciprocal  total  pulse  length 
consisting  of  many  hundreds  of  micro-pulses.  In  this 
paper  we  briefly  review  previous  work  on  the  target 
range-Doppler  image  features  which  are  useful  for 
target  discrimination,  and  then  we  examine  the 
detection  of  the  target  under  low  signal  conditions 
given  passive  sensor  or  radar  handover  range  and 
velocity  information  which  contains  errors.  A 
Gaussian-mode  radiometry  analysis  is  presented,  the 
range-Doppler  image  formation  process  is  reviewed, 
and  the  effects  of  target  speckle,  pointing  errors,  and 
atmospheric  turbulence  are  discussed.  A  metric  for 
determining  a  "target  detection,”  under  these  low- 
signal  initial  conditions,  is  proposed. 

1.0  Introduction 

Under  the  Discriminating  Interceptor  Technology 
Program  (DITP)  and  Advanced  Sensors  Technology 
Program  (ASTP)  at  BMDO^’‘,  ultra-compact  mode- 
locked  CO2  and  mode-locked  Nd:fiber  ladars  are 
being  developed  which,  when  used  in  conjunction 


with  passive  sensors  or  radar,  will  enhance  the 
discrimination  robustness  of  future  interceptors. 
Missile  seeker  technology  development  is  the 
emphasis  of  DITP,  and  surveillance  sensor 
technology  is  the  focus  of  ASTP.  Combining  a  ladar 
and  passive  sensor  outputs  will  increase  the  number 
of  target  features  observed  due  to  the  complementary 
nature  of  the  target  signatures  measured.  For 
example,  ladar  measurements  of  range-resolved 
intensity  and  velocity  are  not  directly  obtainable  with 
passive  sensor  measurements.  Likewise,  in 
conjunction  with  a  radar,  the  ladar  provides  precision 
target  angular  information  as  well  as  more  precise 
range,  Doppler,  and  range-resolved  Doppler 
measurements. 

The  mode-locked  lasers  under  development^’^ 
produce  a  train  of  coherent  narrow  pulses,  or  "micro- 
pulses,"  of  pulsewidths  from  1.3  ns  to  2.9  ns.  This 
corresponds  to  a  down-range  resolution  of  from  20 
cm  (NMD)  to  44  cm  (TMD).  Also,  the  frequency  of 
the  CO2  laser  is  28,300  GHz  (10.6  /im)  for  exo- 
atmospheric  operation  or  26,900  GHz  (11.15  fim 
isotopic  CO2)  for  endo-atmospheric  operation.  The 
frequency  of  the  Nd:fiber  laser  is  283,000  GHz  (1.06 
/im).  Since  target  velocity  resolution  is 
approximately  K/2T^aveform»  we  see  that  the  velocity 
resolution  is  about  3,000  or  30,000  times  higher  than 
that  of  X-band  radar.  (Or,  the  same  measurement 
can  be  done  in  1/3, 000th  to  l/30,000th  of  the  time.) 
The  target  image  Doppler  resolution  is  also 
determined  by  the  reciprocal  total  "macro-pulse"  time 
consisting  of  many  hundreds  of  micro-pulses. 
Consequently,  it  is  possible  to  form  a  precision  target 
range-resolved  Doppler  image  in  real  time  which  is 
useful  for  target  discrimination.  In  this  paper  we  will 
concentrate  on  the  CO2  ladar  systems. 
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2.0  Image  Feature  Extraction  for  Discrimination 

Previous  papers^'^  have  examined  feature 
extraction  from  ladar  range-resolved  Doppler  images 
for  target  discrimination.  Two  range-Doppler 
images'*  of  cones,  taken  at  the  U.S.  Army  AMOR 
facility,  are  shown  in  Figure  1  for  cone  aspect  angles 
of  island  30°.  Target  surface  produced  laser  light 
speckle,  due  to  the  rough  surface  with  respect  to  X, 
can  be  seen  in  the  images.  This  speckle  effect  can 
be  reduced  by  a  non-coherent  summation  of 
consecutive  range-Doppler  images,  although  it  is 
difficult  to  maintain  perfect  range  and  frequency 
registration  for  precise  summing  and  it  requires 
additional  time.  The  target  features  that  can  be 
estimated’-'’  from  the  range-Doppler  image  include: 
target  length,  diameter,  spin  rate,  precession  rate, 
precession  half-angle,  center  of  mass  location,  total 
cross-section,  range-resolved  cross-section,  cross- 
section  fluctuation  statistics  in  time  and  frequency 
domains,  and  aspect  angle  variations  of  total  and 
resolved  cross-section.  A  sub-set  of  these  estimated 
parameters  has  been  found'*  to  be  adequate  for  RV 
and  decoy  discrimination. 

3.0  Gaussian  Mode  Coherent  Ladar  Analysis 


The  ladar  carrier-to-noise  ratio  (CNR),  or  signal- 
power  to  local-oscillator-noise-power  ratio,  is 
determined  by  integrating  the  transmitted  laser  beam 
and  the  back-propagated-local-oscillator  laser  beam 
(including  turbulence  effects  on  both  beams)  over  the 
surface  of  the  unresolved  target.  This  treatment 
follows  the  work  of  Frehlich*,  also  known  as  the 
"beam-overlap"  technique.  We  propagate  the  laser 
beam  from  the  laser  to  the  transmitting  aperture  and 
then  out  to  the  target,  and  we  also  back-propagate  the 
LO  beam  incident  on  the  detector  (including 
clipping),  back  out  to  the  receiving  aperture  (more 
clipping),  and  then  back  out  to  the  target.  This  is 
denoted  as  the  BPLO  mode.  The  two  modes  are  then 
integrated  over  the  target  surface  where  a  diffuse 
target  scattering  function  can  be  used  to  simplify  the 
complicated  integrals.  This  simplification  results  in 
a  CNR  given  by 
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where  It  and  IgpLo  are  the  laser  irradiance  at  the 


transmitter  aperture  propagated  out  to  the  target  and 
the  back-propagated-local-oscillator  beam  irradiance 
on  the  target,  excluding  the  atmospheric  extinction. 
The  brackets  indicate  averaging  over  atmospheric 
turbulence  effects.  "B"  is  the  electrical  bandwidth 
assumed  to  be  flat-top  and  subtending  all  of  the 
received  signal  frequencies  and  yet  no  more  than 
necessary  of  the  local  oscillator  noise  caused  by  Plo 
incident  on  the  detector,  "p"  is  the  unresolved  target 
scattering  function  in  m^sr  described  below,  and  "ti" 
is  the  detector  quantum  efficiency,  0.5  for  COj  and 
0.85  for  Ndrfiber.  "a"  is  the  atmospheric  extinction 
coefficient  often  taken  to  be  that  of  a  mid-latitude 
summer  atmosphere  with  rural  aerosol  with  23  km 
visibility  at  low  altitudes  and  background-plus-aged- 
volcanic  stratospheric  aerosols.  The  integral  over  a 
is  computed  in  spherical-shell  layers  along  the  path 
of  the  laser  mode  and  the  backscattered  speckle- 
pattern  return.  (Note,  a  is  zeroed  for  exo- 
atmospheric  operation.)  Higher-order  terms^  in  the 
integrand  of  (1)  result  in  an  occasional  +3dB 
"backscatter  enhancement"  which  is  neglected  here. 

The  transmitted  and  BPLO  mode  irradiances  at  a 
distance  R  from  the  ladar  are  both  Gaussian: 

'  '  Jt  (o\R) 


in  good  agreement  with  Technology  Development 
Corp.,  Hughes,  and  Textron  Systems  Div.  mode- 
locked  laser  measurements.  We  assume  optical  train 
transmissions,  t,  and  t^,  of  0.7  each.  The  laser  mode 
l/e“  irradiance  radius  at  range  R  is  given  by 

\2 
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which  includes  the  effect  of  atmospheric  turbulence 
on  the  average  irradiance*  through  Po(R).  is  the 
1/e^  irradiance  beam  waist  radius  located  near  the 
aperture  of  the  telescope.  The  optical  train  and 
boundary  layer  aberration  effects  produced  Strehl 
ratio.  Sop,,  is  assumed  to  be  negligible,  equal  to  1.0 
here.  The  transverse  coherence  distance  p„(R)  given 


by*: 


Po(R)  = 


(4) 


where  z  =  0  corresponds  to  the  ladar  transmitting  / 
receiving  aperture.  The  quantity  C„\h)  is  the 
"refractive  index  structure  coefficient"  vs  altitude. 
The  Huffnagle-Valley  5/7  and  clear  one  night  models 
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Figure  1,  AMOR  range-Doppler  image'*  of  a  spinning  cone  at  15°and  30®. 
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igure  2,  Cross-section  (m^/sr)  of  a  10  m  x  1  m  cylinder  with  2  m  conical  tip,  5%  diffuse  reflectivity. 
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Figure  3.  Cross-section  (m7sr)  of  a  2  m  x  1  m  cone,  10%  and  5%  diffuse  reflectivity. 
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are  often  used.  (For  exo-atmospheric  ladars,  =  0 
and  is  infinite.)  We  see  that  the  average  effect  of 
atmospheric  turbulence  is  to  spread  the  transmitted 
and  BPLO  modes  and  thereby  decrease  the  received 
signal  strength.  At  CO2  laser  wavelengths  this  effect 
is  usually  only  a  few  percent  unless  the  TMD  target 
is  still  close  to  the  ground. 

3.1  Gaussian  Mode  Strehl  Ratios 

The  laser  mode  Strehl  ratio,  S^,  is  reduced  by 
telescope  clipping  and  central  obscuration.  Degnan 
and  Klein^  have  shown  that  for  a  given  aperture 
diameter,  d^^^  =  25  cm  here,  the  Strehl  ratio  is 

=  (exp(-a-Y^)-exp(-a2))7a^  (5) 

where 

a  =  clear  aperture  radius  /  (=  1.12  for 

maximum  irradiance  @7  =  0) 

Y  =  linear  obscuration  ratio 

(0^  =  12.5  cm  /  1.12  =  11.125  cm  for  maximum 
irradiance 

full  angle  =  4X/(Jt  2coJ  =  63.8 //r@  11.15 
@  1/e^  irradiance 

Yura®  gave  identical  results  for  y  =  0.  However,  to 
get  the  correct  ’’Strehl  ratio”  for  use  in  equation  (2), 
we  would  ordinarily  solve  equation  (5)  for  the 
optimum  a,  and  then  substitute  this  value  of  a  into 
the  Strehl  ratio  equation  for  a  given  (which  is 
the  correct  S^J: 

L,  "  -exp( -ar)y 

This  is  because  the  far- field  irradiance  distribution  in 
the  absence  of  aperture  clipping  and  obscuration,  in 
the  definition  of  a  Strehl  ratio,  is  a  function  of  (o^ 
only.  For  maximum  on  axis  irradiance,  8,^^  = 
Sbplo  =  from  equation  (6)  for  CO2.  At  1.06  /nn 
and  a  50  jlit  full  angle  divergence,  these  tcnns  are 
almost  1. 

3.2  Pointing  Bias  and  Jitter  Effects 

Yura^  has  developed  an  approximation  for  the 
average  decrease  in  received  irradiance  due  to 
pointing  jitter  and  bias  errors  for  direct  detection. 
For  coherent  detection  the  mixing  efficiency  is  also 
degraded  by  an  off-axis  target  source,  when  there  is 
a  pointing  error,  due  to  both  a  decrease  in  transmitted 
irradiance  and  BPLO  irradiance,  as  indicated  in 
equations  (1)  and  (2)  above.  Tliereforc,  we  define  a 
Gaussian  transmitted  intensity  times  an  off-axis 
mixing  efficiency  degradation  due  to  a  decrease  in 
the  BPLO  irradiance  as 


/(e)  = /^exp(-0^/20^)exp(-072ef) 

and  repeat  Yura’s  calculation  to  find  the  irradiance 
probability  density  function 


-a-b^n 


( 


\^°l 


aV2-l 


/>^6(ln(///J]^g^ 


where  a  and  b  are  defined  w.r.t.  the  Gaussian  mode 

0  ^  ^  ^  ^  _^(9) 


2jto)„  2  a.  e„ 

Normalizing  I  by  to  define  I„,  we  also  have 

The  average  coherent  detection  intensity  and  mixing- 
efficiency -product,  <I„>,  can  then  be  computed  for 
various  Oj  and  0p,: 
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'coh 


a^+2 


expj 


/  \ 

-b  a 


a^+2 


</> 


(11) 


where  Yura’s  "a"  parameter  approximation  for 
coherent  detection. 

The  ladar  fast  steering  mirror  and  slow  steering 
mirror  will  have  bias  and  jitter  errors.  In  addition, 
the  commanded  azimuth  and  elevation  pointing 
angles  will  also  have  bias  and  Jitter  errors.  We 
assume  3  fir  of  bias,  10  fir  of  jitter,  and  co^,  = 
11.125  cm,  and  we  find  =  0.55  in  this  analysis. 


3.3  Target  Cross-section  Models 

For  the  TMD  target  we  will  assume  that  the 
target  is  a  cylinder  10  m  long  and  1  m  in  diameter 
with  a  2  m  long  conical  tip  having  a  diffuse 
reflectivity  coefficient  of  5%  at  11.15  /rm.  In  Figure 
2  is  the  cylinder+cone  cross-section^“  in  m7sr  vs 
aspect  angle.  The  retro-reflective  and  glint 
components  are  omitted  because  they  are  generally 
negligible,  with  respect  to  the  diffuse  component, 
except  near  90°  aspect.  We  also  choose  this  cylinder 
to  be  at  45°  aspect  angle,  thus,  it  is  7  m  in  depth. 
Each  TMD  ladar  micro-pulse  is  2.9  ns  or  87  cm  in 
length.  When  50%  of  the  target  surface  scattering 
points  within  the  micro-pulse  are  replaced,  we  have 
a  new  speckle  source.  Thus,  every  87  cm/2  produces 
a  new  speckle  source,  and  the  missile  cross-section 
may  be  divided  among  about  14  nearly  independent 
target  down-range  coherent  return  cells.  The  cross- 
section  per  independent  range  cell  is  then  taken  as 
“  0.074  (w  Vsr)  /  14  s  p,^  /jc  (12) 
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where  the  (pA/jt)  indicates  an  area  x  reflectivity  /  sr 
for  clarity.  To  detect  a  target,  only  one  or  more  of 
the  down-range  cells  must  exceed  a  threshold.  (At 
larger  aspect  angles  near  90°,  the  total  cross-section 
goes  up,  and  the  number  of  range  coherence  cells 
approaches  one,  greatly  increasing  the  CNR  per  range 
cell.) 

The  NMD  target  is  taken  to  be  a  2  m  long,  1  m 
diameter  cone  with  5%  diffose  reflectivity  at  10.6 
/im  (and  10%  diffuse  reflectivity  at  1.06  jum)  as 
shown  in  Figure  3.  The  micro-pulses  are  taken  to 
1.3  ns  long  for  20  cm  range  resolution,  thus  the  cone 
at  45°  produces  roughly  7  independent  range  cell 
returns.  Again,  we  neglect  the  retro-reflective  and 
glint  components. 


3.4  Unresolved  Target  CNR  Model 

We  can  now  evaluate  the  overlap  integral  of 
equation  (1)  assuming  a  delta-function  for  the 
unresolved  target  angular  (transverse)  subtense 

p(f,d)  = (m^sr)  (13) 

Jt 


and  assuming  some  pointing  error  offset  "a"  from  the 
optical  axis  of  the  LADAR.  Assuming  this  delta- 
function  target  model,  we  get  using  (1),  (2),  and  (13) 

-  (14) 


CNR(t)=. 


hvB 


/  N 

P/ 

r  tar  tar 

4X2 

f^Yp 

/  \ 
-4a  “ 

K  ] 

jt;2co^(/?) 

[co2(R)j 

a. 


for  aligned  transmitter  and  BPLO  modes.  We  have 
defined  a  ’’system  efficiency”  as 

T1.  (15) 


where  and  SgpLo  are  from  Section  3.1,  and  0.83  is 
the  best  mixing  efficiency^  for  a  Gaussian  LO  mode 
and  the  Airy  disk  with  7  =  0.  We  have  also  included 
an  interference  fringe-contrast-ratio  degradation  term, 
ti^(R),  to  account  for  the  decrease  in  fringe  contrast, 
versus  round-trip  delay  time,  due  to  phase-noise 
within  the  LO  laser  source.  In  this  analysis  we  set 
this  term  to  1.0  corresponding  to  a  perfect  LO  laser. 


4.0  Classical  Receiver-Operating-Characteristic 


In  heterodyne  detection,  the  detector  output  is 
bandpass  filtered  around  an  intermediate  frequency 
(IF)  produced  by  mixing  the  backscattered  light 
signal  with  the  LO  laser  to  maximize  CNR  following 
equation  (1).  For  a  pure  speckle  target  both  the 


signal  and  LO  noise  are  Gaussian  random  processes. 
Following  envelope  detection”  theory,  the  rectified 
and  low-pass  filtered  signal  and  LO  noise  voltages 
are  Rayleigh  distributed.  The  probability  of  a  false 
alarm  is  then  the  integral  from  some  threshold 
voltage  to  infinity  over  the  output  voltage  probability 
density  function  (PDF)  in  the  absence  of  signal: 

\dv  =  exp(-v,  /20io) 


P,.  = 


f 

J  exp_ 

^LO  (  ' 


^  -V2 


2  a 


'LO 


Likewise,  the  probability  of  detection  is  the  integral 
from  the  same  threshold  voltage  over  the  output  PDF 
in  the  presence  of  both  noise  and  signal: 


=  exp(-v//2(aio  +a,))  = 


exp(-v,/2a2o(CJVR+l))  = 


_  pccyfi.i)-' 


(17) 


The  same  result  occurs  for  square-law  device 
detertion. 


4.1  Independent  Range  Cell  Detections 

Because  there  are  14  nearly  independent  coherent 
source  cells  on  the  TMD  diffuse  cylinder  target  at  the 
assumed  45°  aspect,  we  have  14  nearly  independent 
measurements,  and  the  probability  of  one  or  more 
threshold  crossing  detections  out  of  n  is  given  by: 


P^(a/n  ofn)  = 

x-m 


xl 


_«! _ 


^  D 


(1  -P^)"**(18) 


with  m  =  1  and  n  =  14.  From  equation  (18)  we  find 
that  we  need  a  per  range  cell  P^  of  0.15  to  achieve 
an  0.9  Pp  for  all  the  14  range  cells.  From  equation 
(17)  we  also  find  that  we  require  a  CNR  of  about  4 
per  source  region  on  the  target. 

For  the  NMD  cone  target  of  roughly  7  cells,  from 
equation  (18)  we  find  we  need  a  per  range  cell  P^  of 
0.28  to  achieve  an  0.9  Pp  for  7  cells.  From  equation 
(17)  we  find  we  require  a  CNR  of  about  6  per  cell  at 
a  lO  "*  PpA-  We  note  that  this  is  an  improvement  with 
respect  to  umesolved  target  detection  (1  range  bin). 


5.0  Speckle  Modulation  Effects 


The  diffuse  target  surface  roughness  produces 
constructive  and  destructive  interference  effects 
known  as  laser  speckle.  The  dominant  speckle 
modulation  effects  are  reviewed  here. 

(1)  Transverse  speckle  lobe  motion:  The 
maximum  velocity  spread  from  a  rotating  cylinder  is 
It  j  -  sin(0)  (19) 

total  spread  ^ 

and,  therefore,  the  total  Doppler  spread  from  a 
rotating  cylinder  is  given  by 
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=  2a)Dsin(e)/X 


(20) 


The  average  speckle  size  back  at  the 
transmitter/receiver  is 


(  spcckui}  ■”  ^ 

i^sp^ctuz)  =  ^sjn(e) 


0.6m 


(21) 


at  11.15  //m  for  a  nominal  1  m  diameter,  10  m  long 
cylinder  at  45°.  The  speckle  velocity  perpendicular 
to  the  axis  of  rotation  is  given  by 


speckle  1 


=  2co/?sin(9) 


(22) 


and  the  transit  time  over  the  aperture  is 
approximately 


T  _  ^speckU^  _  X  _  1 

'  XIT  "  2i>«.sin(6)  ■  aC 

from  equation  (21)  for  nominal  targets  of  interest. 

(2)  The  rotation  of  new  random  phase  scatter 
points  on  the  cylinder’s  surface  into  the  ladar  FOV 
results  in  a  new  speckle  realization.  This  determines 
a  target  coherence  time,  due  to  "speckle 

exchange.”  George^^  has  shown  that  the  cylinder 
scatter-point-exchange  coherence  time  is  given  by 


'^exchg 


(JO  sin(6) 


32  ms 


(24) 


which  is  considerably  shorter  than  a  jt/2  rotation 
period.  We  also  assume  here  that  pitch  and  yaw 
rates  are  much  smaller  than  the  roll  rotation  rate. 


6,0  Range-Doppler  Image  Signal  Processing 

We  assume  for  this  analysis  that  the  local- 
oscillator-laser  is  set  in  frequency  to  bring  the 
received  Doppler  shifted  laser  light  down  to  an 
intermediate  frequency  (IF)  of  roughly  500  MHz  with 
a  bandwidth  of  ±162.5  MHz  at  the  detector  output 
(single-phase).  This  spectrum  is  shown  in  Figure  4 
where  the  2.9  ns  micro-pulse  bandwidth  is  325  MHz 
(250  MHz  @  -3dB)  and  the  spectral  peaks  are 
separated  by  1/40  ns  or  25  MHz.  The  signal  spectral 
peaks  are  above  that  of  a  single  micro-pulse 
spectrum  where  is  the  number  of  micro-pulses  in 
the  train  =  Tj^acro  /  This  computer  run  assumes 

a  15  //s  macro-pulse,  thus,  the  spectral  peaks  are  375“ 
or  141,000  above  that  of  a  single  micro-pulse,  and 
the  per  micro-pulse  CNR  is  lO  "^  or  -40dB.  Note  that 


the  negative  frequency  components  are  the  complex 
conjugate  of  the  positive  frequency  components  for 
a  real  (single-phase)  signal  and  are  not  shown.  The 
detector  bandwidth  is  seen  to  be  1  GHz  from  the 
local-oscillator  noise  spectrum. 

An  electrical-LO  frequency-locked-loop  then 
brings  the  500  MHz  IF  down  to  200  MHz  in  both 
inphase  and  quadrature  (I  &  Q),  shown  in  Figure  5. 
The  bandwidth  of  this  mixer  stage  is  assumed  to  be 
500  MHz  resulting  in  very  little  phase  distortion  of 
the  signal  due  to  the  flat  transfer  function  selected,  a 
5th-order  Butterworth.  At  this  lower  center 
frequency,  the  signal  can  be  A/D  converted  at  a  2 
GS/s  sampling  rate  on  both  I  &  Q. 

6.1  Digital  Matched  Filter 

Next,  a  digital  matched  filter  (I  &  Q)  is 
convolved  with  each  40  ns  (80  points  @  2  GS/s) 
segment  in  the  frequency  domain  using  FFT’s^^’^"^. 
This  digital  matched  filter  is  constructed  by  sampling 
the  transmitted  micro-pulses  and  constructing  an 
"analytic-function”  reference  waveform  also  at  the 
200  MHz  center  frequency: 

(25) 

=  IFT{[1  2 2...2  2 1  0  0...0  0] 

where  the  length  of  the  O’s  and  2’s  is  one-half  the 
length  of  the  stored  waveform  less  1.  If  we  have  a 
perfect  laser,  the  analytic  reference  may  be  stored  in 
the  computer  memory.  This  nearly  perfect  analytic- 
function  digital  filter  produces  the  spectrum  shown  in 
Figure  6.  We  immediately  see  that  the  negative 
frequency  noise  components  are  greatly  reduced. 
We  thus  pick  up  a  factor  of  2  (+3dB)  in  CNR  from 
that  calculated  in  equation  (14)  which  included  a  2  in 
the  denominator  for  positive  and  negative 
frequencies.  The  digital  matched  filter  is  perfectly 
matched  to  the  spectrum  of  the  return  from  a  point 
source.  Monte  Carlo  simulations  from  spinning 
cones  and  cylinders  indicate  that  the  improvement  of 
the  analytic-function  digital  matched  filter  is  actually 
about  a  factor  of  4  (+6dB)  better  than  the  top  hat 
spectral  models  used  in  the  CNR  analytical  equations. 
This  is  particularly  true  for  asymmetric  micro-pulses 
(due  to  higher-order  mode-medium  effects)  in 
producing  a  near  perfect  spectral  filter  w.r.t.  a  top  hat 
analytical  filter. 

6.2  Range-Doppler  Image  Construction 

The  outputs  of  the  matched  filter  are  then 
rearranged  in  a  waterfall  plot  in  sequence,  each 
segment  being  40  ns  long  (Figure  7).  Zeros  are 
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Figure  4.  Spectrum  of  the  detector  output  centered  at  500  MHz.  CNR  =  lO"^  per  micro-pulse. 
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Figure  5.  Down  shifted  (I&Q)  spectrum  centered  at  200  MHz  after  A/D  conversion. 
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Figure  6.  Spectrum  of  the  analytic  digital  matched  filter  output.  Note  tne  negative  frequency  noise 
components  are  greatly  reduced. 
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added  above  and  below  the  waterfall  sequence  prior 
to  FFTing  the  range  bin  columns.  The  number  of 
zeros  is  determined  so  that  the  waveform  resolution 
(1/Tjnacro)  subtends  at  least  4  or  5  FFT  points  so  that 
triangular  or  spiked  peaks  are  avoided,  and  we  obtain 
smooth  peaks.  For  a  40  //s  waveform,  for  example, 
this  would  require  zero  padding  to  4096: 

^  J  (26) 

25M/fz/4096 


as  opposed  to  zero  padding  to  2048.  FFTing  down 
the  columns  and  taking  the  magnitude  squared  is 
proportional  to  finding  the  power  spectral  density  of 
this  range  bin.  From  spectral  analysis  theory^^  we 
find  the  mean  square  voltage  down  a  range  bin 
column  to  be 

(v2>  =  (27) 

where  the  v-  are  the  element  voltages,  is  the 
sample  time,  is  the  data  length  (N^^),  and 
T 

= _ L.FT(v)FT(vy  (28) 


is  the  "periodogram"  estimate  of  a  power  spectral 
density  (PSD)  with  N,^  equal  to  plus  zero 
padding.  Consequently,  the  average  noise  PSD  in  a 
range-Doppler  image  column  is 
(|v^p>  40ns 


>  = 


4096 


IL  (v^Hz)  (29) 


The  average  signal  PSD  is 


>  = 


40n5  N 


MP 


(1/7’_J/(25M//z/4096) 


,  (30) 

(vVffi) 


when  the  target  bandwidth  is  less  than  the 
waveform  resolution,  as  the  signal  is  only  in  a  few 
Doppler  bins  (the  denominator  of  (30)).  Taking  the 
ratio  (30)  to  (29),  we  find  the  range-Doppler  image 
column  CNR  to  be 


40a75 


CNR^f  (31) 


when  or  <  1  /  One  such 

range-Doppler  image  of  a  2  m  cylinder  at  45°  is 
shown  in  Figure  8  at  a  CNR  of  5.3  after  processing 
375  micro-pulses  (15  jus). 

When  the  target  Doppler  spread  exceeds  the 
waveform  frequency  resolution,  the  target  energy 
gets  spread  into  additional  cross-range  FFT  points. 
For  this  situation  the  average  signal  PSD  in  those 
bins  containing  signal  becomes 


(0>  = 


(|v  p)40nsA/,^ 
6v,.,,«/(25M//z/4096) 


and  taking  the  ratio  of  (32)  to  (29) 

CNR 

= 


MF 


40  ms  6v 


target 


(32) 


(33) 


when  or  >  1  /  Thus,  the 

CNR  remains  constant  once  T^^acro  =  1  /  ^^target  is 
reached,  since  is  fixed.  For  longer  macro¬ 

pulse  times,  non-coherent  effects  come  into  play. 


6.3  Non-coherent  Processing  of  the  Range- 
Resolved  Spectra 

To  further  increase  the  image  CNR,  we  might 
non-coherently  sum  the  range-bin  PSDs  calculated 
every  T^^^  ^  pulse  train.  In  actuality,  the  PSD 

estimator  given  in  equation  (28)  automatically  does 
this.  For  example,  under  the  AW  ACS/EAGLE 
program^^  it  was  determined  that  the  TMD  targets  of 
interest  usually  rotate  at  less  than  1/20  revolution  per 
second.  This  produces  a  Doppler  spread  of  less  than 
2  X  28  kHz  (CO2  wavelengths)  at  90°  aspect  and  less 
than  2  X  19  kHz  at  the  assumed  45°  aspect.  Due  to 
cylindrical  surface  weighting,  the  -3dB  Doppler 
spread  is  around  19  kHz  total,  hence,  we  are  limited 
to  macro-pulse  lengths  of  roughly  54  f4s  (1350  micro¬ 
pulses)  for  efficient  coherent  signal  processing  where 
equation  (31)  applies.  If  we  select  a  long  pulse  train 
of  T^3^^  =  10  ms,  we  then  have  10  ms  /  54  //s  =  185 
independent  spectral  measurements.  If  we  then  non- 
coherently  summed  these  consecutive  independent 
spectra,  we  have,  from  classical  microwave  radar 
analysis^^ 


CNR^,  =  CNR^r 


/  \ 
T  , 

^  j- 

coh 

macro 

(^40/15^ 

T  . 

coh  y 

(34) 


where  the  exponent  actually  varies  from  0.6  to  0.8. 
In  the  formation  of  the  PSD  via  equation  (28)  down 
each  range-bin,  that  is  what  happens.  This  amounts 
to  a  factor  of  up  to  52,000  for  a  10  ms  train  and  a  54 
/IS  target  coherence  time.  It  does,  however,  require 
good  range-bin  and  good  Doppler-bin  tracking  (or 
registration)  in  order  to  produce  an  accurate, 
unblurred  image.  We  note  that  retargeting  rates  of 
50  objects  per  second  would  allow  up  to  a  20  ms 
macro-pulse  length,  further  enhancing  the  summation. 

Assuming  a  10  ms  macro-pulse  train  of  7.4  /J 
micro-pulses  (185  W  average),  a  54  //s  target 
coherence  time,  and  a  +6dB  analytic  digital  matched 
filter  enhancement,  the  points  where  the  TMD 
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Figure  7.  Waterfall  plot^'*  of  the  sequence  of  micro-pulses  after  matched  filtering.  Note  the  slowly 
varying  small  signals  at  the  cone  tip,  and  the  faster  varying  strong  signals  near  the  base. 
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Figure  8.  Low  CNR  range-Doppler  image  of  a  cylinder  2  m  long  at  45°  with  ±26  kHz  Doppler  spread. 
Contour,  Doppler  axis,  range  axis,  and  isometric  views.  CNR  =  3.6x10*^  per  micropulse,  5.3  after 
processing  15  /is  (375  micro-pulses)  of  data.  A  threshold  crossing  of  the  range-axis  view  is  a  detection. 
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CNRrds  reaches  4  are  plotted  in  Figure  9  for 
various  cylinder+cone  target  altitudes.  At  the 
baseline  target  altitude  of  20  km,  we  see  that  the  per 
range-cell  CNR  of  4  is  met  at  550  km  (indicating  a 
full  target  (14  cells)  Pp  of  0.9  at  a  Pp^  =  10  ''  per 
Section  4.  Atmospheric  extinction  at  low  altitudes 
dominates  the  losses.  For  the  exo-atmospheric  NMD 
scenario,  X.  =  10.6  /rm  at  4  /tJ  per  pulse  and  we  see 
the  1/R‘'  dependence  in  CNR^pj  in  Figure  10. 
Maximum  range  (CNR^p^  =  occurs  at  about  370 
km  with  the  cone  target  at  45°  aspect. 

6.4  Alternative  Target  Detection  Algorithms 

Summing  down  the  range  columns  of  the  range- 
Doppler  image  to  form  an  "amplitude  (or  amplitude^) 
vs  range"  image  was  found  to  be  not  as  effective  an 
algorithm  as  a  range-axis  view  of  the  range-Doppler 
image  (magnitude)^.  This  is  because  the  cross-range 
cells  outside  the  target  contain  only  LO  noise  and 
corrupt  the  cross-range  sums.  The  cross-range  bin 
number  must  be  maintained  fairly  large  to 
compensate  for  frequency  errors  in  estimating  the 
target’s  center  of  mass  Doppler  and  in  the  frequency- 
locked  loop.  Likewise,  summing  across  the  Doppler 
bins  to  form  a  "Doppler  resolved  target"  image  is  not 
as  good  for  detection. 

7.0  Work  in  Progress 

1)  The  usual  slowly  spinning  targets  do  not 
generate  Gaussian,  Rayleigh,  negative-exponential 
statistics  when  T^^^  <  T^^-  Monte  Carlo  runs 
indicate  Rician  and  modified-Rician  (sinewave  plus 
diffuse)  statistics  are  observed  with  the  same  <V3“> 
target  signal  magnitude  but  with  a  reduced  variance. 
This  means  the  detection  statistics  will  be  better. 
The  allocation  of  the  sinewave  component  to  the 
diffuse  component  must  be  determined  as  a  function 
of  target  spin  rate  such  that  the  total  target  cross- 
section  sum  remains  equal  to  the  target  true  cross- 
section.  The  detection  of  this  signal  in  Gaussian  LO 
noise  must  then  be  computed.  When  T^^^^  is  greater 
than  T^i,»  Gaussian,  Rayleigh,  negative-exponential 
statistics  return. 

3)  A  target  aspect  angle  closer  to  90°  would 
probably  be  more  appropriate  giving  larger  signal  per 
range  cell. 

4)  Using  the  full  macro-pulse  waveform  as  a 
digital  matched  filter  and  then  forming  a  range- 
Doppler  image  (with  no  Doppler  cross-range 
information)  produces  a  superior  detection 
measurement  compared  to  the  range-axis  view  of  the 


range-Doppler  image.  The  operations  per  second 
required  for  these  long  vectors  may  be  prohibitive. 

5)  When  there  is  a  range  (return  time)  error,  the 
target  image  shifts  forward  and  backward  in  the 
range-Doppler  matrix  not  causing  much  of  a  problem. 
When  there  is  a  center-of-mass  frequency  error,  the 
peak  of  the  matched  filter  output  decreases  and 
distorts,  an  effect  known  as  "range-Doppler 
coupling."  The  magnitude  of  this  effect  is  under 
investigation. 

6)  Sinsusoidal  and  linear  frequency  modulation 
waveforms  are  being  investigated  for  the  endo- 
atmospheric  (ASTP)  applications  as  well. 
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